Massive Black Holes in Star Clusters. II. Realistic Cluster Models 
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ABSTRACT 

We have followed the evolution of multi-mass star clusters containing massive central black 
holes through coUisional iV-body simulations done on GRAPE6. Each cluster is composed of 
between 16,384 to 131,072 stars together with a black hole with an initial mass of Mbh — 
IOOOMq. We follow the evolution of the clusters under the combined influence of two-body 
relaxation, stellar mass-loss and tidal disruption of stars by the massive central black hole. 

We find that the (3D) mass density profile follows a power-law distribution p ^ r^" with slope 
a = 1.55 inside the sphere of influence of the central black hole. This leads to a constant density 
profile of bright stars in projection, which makes it highly unlikely that core collapse clusters 
contain intermediate-mass black holes. Instead globular clusters containing massive central black 
holes can be fitted with standard King profiles. Due to energy generation in the cusp, star clusters 
with intermediate-mass black holes (IMBHs) expand. The cluster expansion is so strong that 
clusters which start very concentrated can end up among the least dense clusters. The amount of 
mass segregation in the core is also smaller compared to post-collapse clusters without IMBHs. 
Most stellar mass black holes with masses AIbh > 5Mq are lost from the clusters within a few 
Gyrs through mutual encounters in the cusp around the IMBH. 

Black holes in star clusters disrupt mainly main-sequence stars and giants and no neutron 
stars. The disruption rates are too small to form an IMBH out of a Mbh ~ 50 Mbh progenitor 
black hole even if all material from disrupted stars is accreted onto the black hole, unless star 
clusters start with central densities significantly higher than what is seen in present day globular 
clusters. 

We also discuss the possible detection mechanisms for intermediate-mass black holes. Our 
simulations show that kinematical studies can reveal 1000 Mq IMBHs in the closest clusters. 
IMBHs in globular clusters are only weak X-ray sources since the tidal disruption rate of stars 
is low and the star closest to the IMBH is normally another black hole which prevents other 
stars from undergoing stable mass transfer. For globular clusters, dynamical evolution can push 
compact stars near the IMBH to distances small enough that they become detectable sources 
of gravitational radiation. If 10% of all globular clusters contain IMBHs, extragalactic globular 
clusters could be one of the major sources of gravitational wave events for LISA. 

Subject headings: black hole physics — globular clusters — methods: N-body simulations — stellar dynamics 



1. Introduction 

This is the second paper in a series of iV-body 
simulations which deal with the dynamical evo- 
lution of star clusters containing massive central 
black holes. In Baumgardt, Makino & Ebisuzaki 
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(2004) (paper I), we followed the evolution of sin- 
gle mass clusters with black holes that contained a 
few percent of the total system mass. We showed 
that the density distribution of stars inside the 
sphere of influence of the black hole follows a 
p ^ 7.-1-75 power law, in good agreement with 
results from previous analytical estimates and in- 
direct simulation methods (Bahcall & Wolf 1976; 
Cohn & Kulsrud 1978; Marchant & Shapiro 1980). 
We also derived the rate of tidal disruption of 



stars. The present paper extends these resuhs to 
the more reahstic but much less studied case of a 
star cluster with a spectrum of stellar masses and 
is aimed at globular clusters with intermediate- 
mass black holes of a few hundred to a few thou- 
sand solar masses at their centres. 

X-ray observations of starburst and interacting 
galaxies have revealed a class of ultra-luminous X- 
ray sources (ULX) , whose luminosities exceed the 
Eddington luminosities of stellar mass black holes 
by orders of magnitude (Makishima et al. 2000), 
making them good candidates for IMBHs. The ir- 
regular galaxy M82, for example, hosts an ULX 
with luminosity L > 10^° ergs/sec near its cen- 
tre (Matsumoto et al. 2001; Kaaret et al. 2001) 
the position of which coincides with that of the 
young (T « 10 Myrs) star cluster MGG-11. Porte- 
gies Zwart et al. (2004) have performed A^-body 
simulations of several star clusters in M82, using 
the cluster parameters determined by McCrady et 
al. (2003). They found that runaway merging of 
massive stars could have led to the formation of 
an IMBH with a few hundred to a few thousand 
solar masses in MGG-11, thereby explaining the 
presence of the ultraluminous X-ray source. Gen- 
eral conditions when runaway merging of stars can 
lead to the formation of IMBHs were discussed in 
Portegies Zwart & McMillan (2002) and Rasio et 
al. (2004). 

Apart from run- away collisions of main se- 
quence stars, massive black holes could also be 
build up through the merger of stellar mass black 
holes (Mouri & Taniguchi 2002) via gravitational 
radiation in dense enough star clusters. In less 
concentrated clusters, this process is also possi- 
ble but may take up to a Hubble time (Miller & 
Hamilton 2002). Finally, intermediate-mass black 
holes could also form by the accretion onto a 
stellar-mass black hole of interstellar cluster gas 
due to radiation drag from bright stars, provided 
the stellar mass function of the cluster is shallow 
enough (Kawakatu & Umcmura 2004). 

Although there are several possible ways to 
form IMBHs, the observational evidence for 
intermediate-mass black holes in star clusters is 
as yet much less clear. Due to their high cen- 
tral densities and since they are relatively close, 
galactic core-collapse clusters are obvious places 
to look for intermediate-mass black holes. In- 
deed, for almost 30 years, M15 has been thought 



to harbour an intermediate-mass black hole in 
its centre (Newell et al. 1976). The most recent 
analysis of this cluster was done by Gerssen et al. 
(2002, 2003), who used HST to obtain new spec- 
tra of stars in the central cluster region. They 
found that the velocity dispersion can best be 
explained by an intermediate-mass black hole of 
mass (1.7±2.7) • IO'^Mq. However, Baumgardt et 
al. (2003a) have shown that the observational data 
can also be explained by the core-collapse profile 
of a 'standard' star cluster without a massive cen- 
tral black hole. In this case the central rise of the 
mass-to-light ratio is created by an unseen con- 
centration of neutron stars and heavy mass white 
dwarfs. Such a model is also able to explain the 
velocity dispersion derived from the proper motion 
of stars near the centre of M15 (McNamara et al. 
2003, 2004). Similarly, a dense concentration of 
compact remnants might also be responsible for 
the high mass-to-light ratio of the central region 
of NGC 6752 seen in pulsar timings (Ferraro et al. 
2003; Colpi et al. 2003). Outside our own galaxy, 
Gebhardt et al. (2002) have reported evidence for 
a 20,000 Mq black hole in the M31 globular clus- 
ter Gl, but Baumgardt et al. (2003b) showed that 
dynamical simulations without black holes com- 
pletely explain the observed velocity dispersion 
and density profile of this cluster. 

Despite the unclear observational situation, the 
presence of intermediate-mass black holes in star 
clusters remains an interesting possibility. They 
would provide the missing link between the stel- 
lar mass black holes which form as a result of 
stellar evolution and the supermassive black holes 
in galactic centres (Ebisuzaki et al. 2001). They 
would also be prime targets for the forthcoming 
generation of both ground and space-based grav- 
itational wave detectors due to their high masses 
and the fact that they have the potential to merge 
with other black holes if residing in dense star clus- 
ters. 

In this paper we study the dynamical behaviour 
of massive black holes in globular clusters. Section 

2 describes the set-up of our runs and in section 

3 we present our main results concerning the dy- 
namical evolution of the star clusters, the tidal 
disruption of stars and the possibilities of detect- 
ing the central black hole. In section 4 we report 
our conclusions. 



2. Description of the runs 

We simulated the evolution of clusters contain- 
ing between A^ = 16, 384 and 131, 072 (128K) stars 
using Aarseth's coUisional A^-body code NB0DY4 
(Aarseth 1999) on the GRAPE6 computers of 
Tokyo University (Makino et al. 2003). Most clus- 
ters were treated as isolated and followed King 
Wq ~ 7.0 profiles initially. Simulations were 
run for a Hubble time, which was assumed to be 
T = 12 Gyrs. Since black holes probably form 
early on in the evolution of a globular cluster by 
e.g. run-away merging of massive stars, we started 
our runs with a massive black hole at rest at the 
cluster centre. As in paper I, we modified the ve- 
locities of the cluster stars to prevent the cluster 
centre from collapsing after adding the IMBH to 
the cluster. 

All clusters started with a central black hole 
of mass Mbh = IOOOMq. If the Mq - a relation 
found by Gebhardt et al. (2000) for galactic bulges 
holds for globular clusters as well, this corresponds 
to the IMBH mass expected in a typical globular 
cluster. In addition, the best case found so far 
for an IMBH in a star cluster, M82 X-1 in MGG- 
11, must have a mass between a few hundred to 
a few thousand solar masses based on its X-ray 
luminosity and the frequency of quasi-periodic os- 
cillations seen in the X-ray flux (Matsumoto et al. 
2001; Strohmayer & Mushotzky 2003). 

Stellar evolution was treated by the fitting for- 
mulae of Hurley et al. (2000), assuming a mean 
cluster metallicity of Z = 0.001 and a neutron star 
retention fraction of 15%. The 15% retention frac- 
tion was imposed by immediately removing 85% 
of all neutron stars at the time of their forma- 
tion while leaving the velocities of the remaining 
ones unchanged. The exact form of the IMF at 
the low-mass end will not critically influence the 
results of our simulations as such stars are mere 
test particles in the gravitational field of the higher 
mass stars. More important is the mass function 
at the high-mass end, especially the number and 
mass distribution of black holes formed during the 
run. Unfortunately, the initial-to-final mass rela- 
tion for high-mass stars is currently not precisely 
known since it depends among other things on the 
assumed amount of stellar wind mass-loss in the 
final phases of stellar evolution and the details of 
the explosion mechanism (Fryer & Kalogera 2001). 



In addition the metallicity of the progenitor star 
will effect the mass of the final black hole (Heger et 
al. 2004) for low-metallicity stars, so globular clus- 
ters with different metallicities might have differ- 
ent black hole mass distributions. Since the frac- 
tion and mass distribution of heavy mass black 
holes could have a strong influence on the out- 
come of our simulations, we performed two series 
of simulations. 

In our first series of simulations, the mass func- 
tion of the cluster stars was given by a Kroupa 
(2001) IMF with a lower mass limit of O.IMq 
and an upper mass limit of 30Mq. Since this 
upper mass limit is only slightly above the mass 
where black holes instead of neutron stars form 
due to stellar evolution, these clusters contain only 
a small fraction of stellar mass black holes, all of 
them with masses below 3Mq. In the second se- 
ries of simulations we used an upper mass limit 
of 100 Mq and transformed the stars directly into 
black holes without further mass-loss. In this se- 
ries, the most heavy mass black holes formed have 
about 45 Mq. In both types of simulations, all 
black holes were retained in the clusters. These 
two cases are the most extreme models, and most 
likely real globular clusters fall in between. 

We did not include a primordial binary popu- 
lation. The presence of a primordial binary popu- 
lation might help the formation of IMBH through 
enhancing the stellar collisions (e.g. Fregeau et 
al. (2003)). Its effect on the structure of a cluster 
with central IMBH would be to decrease the cen- 
tral density through hardening. Thus, our present 
result without primordial binaries gives the upper 
limit of the central density. 

Stars were assumed tidally disrupted if their 
distance to the central black hole was smaller than 
the critical distance given by eq. 3.2 in Kochanek 
(1992): 
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where Mbh and M* are the mass of the black 
hole and the star and i?* is the stellar radius. The 
stellar radii were taken from the formulae of Hur- 
ley et al. (2000). The mass of tidally disrupted 
stars was added to the mass of the central black 
hole. So far we have not incorporated the effects of 
gravitational radiation in our runs since the cen- 
tral densities reached in our simulations are not 



Table 1: Details of the performed A^-body runs. 
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6.23 
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1000.0 
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12 



Notes: 

1. For clusters with more than one simulation, parameters given are average values. 



large enough that gravitational radiation becomes 
important. We also did not include stellar colli- 
sions, which become dynamically important when 
the velocity dispersion around the black hole be- 
comes equal to the escape velocity from the stellar 
surface, as stars cannot undergo large angle en- 
counters any more (Frank & Rees 1976). Even for 
main sequence stars, this corresponds to distances 
of rcoii = 10~^ pc from the black hole, which is 
far inside the distance of the innermost stars from 
the IMBH. 

Table 1 gives an overview of the simulations 
performed. It shows the number of cluster stars 
A'', the number of simulations Nsim performed for 
a given model, the initial concentration Wq of the 
King model, the upper mass limit of the IMF, ini- 
tial cluster mass and half-mass radius. Also shown 
are the final half-mass radius, and the black hole 
masses at the start and the end of the runs. The 
final column gives the (average) number of tidal 
disruptions. 

3. Results 

3.1. Cluster expansion 

Fig. 1 depicts the evolution of lagrangian radii 
for the first four clusters of Table 1. All clusters 
expand since mass-loss of individual stars due to 
stellar evolution decreases the potential energy of 
the clusters and two-body processes in the cusp 
around the black hole exchange energy between 



the stars. As a result, the innermost stars are 
pushed to more negative energies until they are 
tidally disrupted by the IMBH while the rest of the 
cluster stars gain energy and the cluster expands. 
This behaviour is similar to the one seen in the 
single mass runs of paper I (Fig. 5). 

For all clusters the expansion is strongest in the 
initial phase since the cluster radii and therefore 
the two-body relaxation time is smallest in the 
beginning. In addition, the mass-loss of stars due 
stellar evolution is strongest within the first Gyr. 
The expansion is smaller for high-A^ clusters due 
to their longer relaxation times. Table 1 shows 
that for a given mass of the central black hole, the 
final half-mass radius depends only on the num- 
ber of cluster stars and is nearly independent of 
the cluster composition and initial half-mass ra- 
dius. For the N=16K clusters for example, the 
final half-mass radius changes by less than 10% if 
the initial half-mass radius is reduced by a factor 
of 6. The final radius also does not depend much 
on the initial IMF or the initial concentration of 
the King model. 

The two-body relaxation time of a cluster with 
mass Mc and radius rh is given by (Spitzer 1987): 
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for a cluster with stars of average mass m. If 
the cluster expansion is caused by two-body re- 
laxation, half-mass radii of clusters with masses 
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Fig. 1. — Lagrangian radii as a function of time 
for the first four cluster simulations with particle 
numbers between 16,384 < iV < 131,072. All 
clusters expand as a result of energy generation 
in the cusp around the black hole and initially 
also due to stellar evolutional mass-loss which de- 
creases the binding energy of the cluster. Since 
two-body relaxation drives the expansion, high-A'^ 
models expand less than low-A ones. 

— 1/3 

Mc should satisfy a relation r^ ^ M^ after a 
long enough time has passed since the half-mass 
relaxation time is the same for all clusters in this 
case so they expand with the same rate. The fi- 
nal half-mass radii of the clusters in Table 1 are in 
good agreement with such a scaling law. 

Fig. 2 compares the projected half-light radii 
of the clusters in Fig. 1 with the projected half- 
light radii of galactic globular clusters. The re- 
sults of the A^-body simulations are well fitted by 

— 1/3 

an Thp ~ M(j law (solid line), so projected half- 
mass and half-light radii follow the same relation 
as the 3D ones. Projected half-light radii of galac- 
tic globular clusters are taken from Harris (1996). 
Masses for globular clusters were calculated from 
their absolute V magnitudes, assuming a mass-to- 
light ratio of ra/Lv — 2.0. Most globular clusters 
have projected half-light radii which are smaller 
than the ones predicted by an extrapolation of our 
runs. 

Rasio et al. (2004) and Portegies Zwart et al. 
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Fig. 2. — Projected half-light radii of globular 
clusters Thp against their mass (triangles). Pro- 
jected half-light radii of the clusters in the A^- 
body simulation are marked by filled circles. They 
follow a relation rhp ^ M^ (solid line), in 
agreement with the idea that two-body relaxation 
drives the cluster expansion. The expansion due 
to an IMBH is strong enough that clusters with 
IMBHs can end up among the least concentrated 
clusters. 



(2004) have shown that core-collapse times of less 
than a Myr are necessary to form an IMBH in the 
centre of a star cluster by run-away collision of 
main sequence stars. For high-A^ clusters this cor- 
responds to central densities of p,, = IO^Mq/pc"^ 
or higher. Most other processes which have been 
proposed as formation mechanisms for IMBHs also 
require high density environments. Such densi- 
ties are among the highest found in globular clus- 
ters (see Table 2 of Pryor & Meylan 1993). Figs. 
1 and 2 show that even if clusters with massive 
black holes start with very high densities, the sub- 
sequent cluster expansion increases their radii by 
such an amount that they can end up among the 
least dense clusters. Thus, the low density of 
present-day globular clusters do not rule out the 
formation of IMBH, since they might have been 
much more compact when they were born. 

The fact that most clusters have half-light radii 
below our prediction does not speak against the 



presence of lA^EBHs in these clusters, since, while 
our clusters were isolated, galactic globular clus- 
ters are surrounded by a tidal field, which limits 
their growth. Tidally limited clusters which lie be- 
low our predicted line might therefore still contain 
IMBHs. 

The half-mass radius of Gl, the most massive 
globular cluster in M31, is rf^ — 8 pc (Baumgardt 
et al. 2003b). If Gl started much more concen- 
trated and its current size is due to an expansion 
similar to the one seen in our runs, the mass of 
the central IMBH must be larger than 1000 Mq 
since the half-mass radius of Gl is much higher 
than predicted by our runs. The same is true for 
CO Cen, the most massive galactic globular cluster. 

3.2. Density profile 

In paper I it was shown that the density profile 
of a single-mass cluster follows a power-law profile 
p ^ r^" inside the influence radius of the black 
hole with slope a = 1.75, in agreement with results 
from Fokkcr-Planck and Monte Carlo simulations. 
It was also shown that the sphere of influence of 
the black hole is limited by two conditions. For 
small-mass black holes, clusters have a near con- 
stant density core outside the sphere of influence 
of the black hole and the central cusp extends only 
up to a radius r^ at which the velocity dispersion in 
the core becomes comparable to the circular veloc- 
ity of stars around the black hole. The following 
relation was found to give a good estimate for r^: 



GM 



BH 



2 <v2> 



(3) 



where Vc is the core velocity dispersion. For black 
holes which contain more than a few percent of 
the cluster mass, a second condition for ri was 
found to be that the mass in stars inside r^ must 
be smaller than the mass of the central black hole, 
since otherwise the self-gravity of stars becomes 
important and changes the density law. 

Fig. 3 depicts the final density profile of the first 
four clusters after 12 Gyrs. Shown is the three- 
dimensional mass density of all stars. In order 
to calculate the density profile, we have superim- 
posed between 5 (128K) to 20 (16K) snapshots 
centered at T=12 Gyrs, creating roughly the same 
statistical uncertainty for all models. All snap- 
shots were centered on the position of the IMBH. 
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Fig. 3. — 3D mass density profile after T=12 
Gyrs for 4 cluster simulations starting with par- 
ticle numbers between 16,384 < A^ < 131,072. 
Solid lines mark the A^-body results, dashed lines 
a single power-law fit to the density profiles inside 
the radius of influence of the black hole (shown 
by a solid circle). For all models we obtain slopes 
near a = 1.55 for the central stellar cusp. 

We then fitted the combined density profile in- 
side the influence radius of the black hole with 
a power-law density profile. It can be seen that 
we obtain a power-law profile inside r^ with slope 
around a =^ 1.55. There is no visible dependence 
on the particle number. It will be shown in section 
3.3, where the mass segregation in the clusters is 
discussed, that the most massive stars in our clus- 
ters still follow an a = 1.75 power- law, but that 
they are not numerous enough to dominate the 
central region which is the reason for the flatter 
overall slope seen in our runs. 

For N = 16K, the black hole contains more 
than 10% of the total cluster mass at T = 12 Gyrs 
and it dominates the density proflle throughout 
the core. For clusters with particle numbers more 
realistic for globular clusters, the central black 
holes contains an increasingly smaller fraction of 
the total cluster mass, so the velocity criterion lim- 
its the influence of the black hole. In this case the 
central cusp contains only a fraction of the mass 
of the central black hole, for N = 128K for exam- 
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Fig. 4. — Projected density profile of bright stars 
(top) and projected velocity dispersion of the clus- 
ter simulation starting with TV = 131,072 stars. 
The projected distribution of bright stars has a 
constant density core, similar to that seen in most 
globular clusters. Observations of the velocity dis- 
persion could reveal the black hole if a sufRciently 
large number of stars at radii r/rh < 0.01 can be 
observed (bottom panel). 

pie only about 10%, i.e. IOOM0. Since a consid- 
erable fraction of these stars would not be easily 
visibly to an observer because they are compact 
remnants and therefore too faint, the direct ob- 
servation of this cusp for globular clusters with 
Mbh < IOOOMq IMBHs is nearly impossible due 
to statistical uncertainties. 

The upper panel of Fig. 4 depicts the projected 
distribution of bright stars for the cluster with 
N = 128K stars. We define bright stars to be all 
stars with masses larger than 90% of the mass of 
turn-off stars which are still main-sequence stars 
or giants at T = 12 Gyrs. Their density distribu- 
tion should be representative of the distribution 
of light in the cluster. The projected density dis- 
tribution of bright stars does not show a central 
rise and can instead be fitted by a model with a 
constant density core according to 
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{l + r/aY 
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where po and a are constants. A cluster with a 
massive central black hole would therefore appear 
as a standard King profile cluster to an observer, 
making it virtually indistinguishable from a star 
cluster before core collapse. Core collapse clusters 
have central density profiles that can be fitted by 
power-laws with slopes of a ~ 0.7 (Lugger et al. 
1995), which is in contradiction with this profile. 
Since the central relaxation times of core collapse 
clusters are much smaller than a Hubble time, any 
cusp profile would have been transformed into a 
constant density core if an IMBH would be present 
in any of these clusters (see Fig. 5). The pres- 
ence of IMBHs in core collapse clusters is there- 
fore ruled out unless their composition is radically 
different from our clusters. 

The lower panel of Fig. 4 shows the velocity 
dispersions, both the measured one and the one 
inferred from the mass distribution of stars. The 
inferred velocity dispersions were calculated from 
Jeans equation (Binney & Tremainc 1986, eq. 4- 
54) and different mass distributions under the as- 
sumption that the velocity distribution is isotropic 
(i.e. /? = 0). The velocities calculated from the 
mass distribution of the cluster stars alone give 
a good fit at radii r/rh > 0.2 where the mass in 
stars is dominating (except at the largest radii, 
where the velocity distribution becomes radially 
anisotropic). At radii r/r^ < 0.2, the contribu- 
tion of the black hole becomes important. At a 
radius r/r^, = 0.01, the velocity dispersion is al- 
ready twice as high as the one due to the stars 
alone. For a globular cluster at a distance of a few 
kpc, such a radius corresponds to central distances 
of one or two arcseconds. Of order 20 stars would 
have to be observed to detect the central rise at 
this radius with a 95% confidence limit. This 
seems possible both for radial velocity (Gerssen 
et al. 2002) or proper motion studies (McNamara 
et al. 2003) with HST. For a nearby globular clus- 
ter, the detection of a massive central black hole 
through kincmatical studies is therefore possible. 
Similarly, Drukier & Bailyn (2003) concluded that 
the IMBH can be found by studying the tail of the 
velocity distribution through proper motions. 

Fig. 5 shows the projected density distribution 
of stars for the N — 128K cluster at four different 
times. In the beginning, the cluster has a con- 
stant density core due to the initial King model. 
As the cluster evolves, a central cusp forms around 



10" 

106 

§ 10" 

1- 

I 10= 
a 10 

1 

10= 
10' 
S 10" 
10' 
10= 
10 

1 



tt. 



N 



T = Myrs 



♦-• r^*--.^ 



^- T = 200 Myrs 




0,01 0,1 

K/K„ 



0,01 0,1 1 

K/R„ 



Fig. 5. — Projected density distributions of stars 
for four different times for the cluster with A^ = 
131, 072 stars. The densities of bright stars and of 
aU stars are shifted to match each other at i? = 
IORh- The density distribution of aU stars devel- 
ops a cusp profile after a few 100 Myrs. Bright 
stars show a constant density core throughout the 
evolution. 
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Fig. 6. — Average mass of stars in cluster simula- 
tions with different particle numbers. In a cluster 
without an IMBH (dotted line, N=64K) the av- 
erage mass in the core increases until it reaches a 
maximum at core-collapse time (T — 10.5 Gyrs). 
In contrast, the average core mass for clusters with 
IMBHs stays nearly constant throughout the evo- 
lution at <m>— 0.6, independent of N . 



the black hole. This cusp is visible in projection 
already after T ~ 200 Myrs and is fully devel- 
oped around T = 1 Gyrs. Since both times are 
much smaller than a Hubble time, globular clus- 
ters which formed an IMBH early on in their evo- 
lution have reached an equilibrium profile in their 
centre. Throughout the evolution, the distribu- 
tion of bright stars displays a constant density core 
since mass segregation leads to an enhancement of 
heavy-mass compact remnants in the centre. Ex- 
cept for the very first phases after IMBH formation 
when the density distribution might differ from the 
King profiles with which we started, clusters with 
massive black holes should exhibit constant den- 
sity cores in their light profile. 

3.3. Mass segregation 

Due to relaxation, massive stars sink into the 
centre of a star cluster in order to achieve energy 
equipartition between stars of different masses 
(Spitzer 1969). Baumgardt & Makino (2003) and 
Giirkan et al. (2004) showed that mass segregation 



of high-mass stars proceeds on the same timescale 
as the evolution of the clusters toward core col- 
lapse. Baumgardt & Makino (2003) also found 
that by the time a globular cluster goes into core 
collapse, the majority of stars in the core are com- 
pact remnants. 

Fig. 6 shows the average mass of stars in the 
core (defined to contain the innermost 3% of the 
cluster mass) and the average mass of all cluster 
stars for the first 4 clusters of Table 1. The average 
mass of all cluster stars (dashed lines) drops due to 
stellar evolution, which is most effective within the 
first Gyr. For low- TV models the core mass rises 
initially since the relaxation time is short enough 
to allow heavy main sequence stars to spiral into 
the core. For higher- A^ models, the relaxation time 
is larger than the stellar evolution time and stars 
with masses M > 5Mq turn into compact rem- 
nants and lose a large fraction of their mass before 
reaching the core. 

For N=64K stars, we also performed a compari- 
son run which started from a King Wq = 7.0 model 




Fig. 7. — Density of different mass groups as func- 
tion of radius for the N — 128K star cluster at 
T=12 Gyrs. Stars in the most massive group fol- 
low a distribution that drops as a « 1.75 in the 
cusp around the black hole. For lower mass stars, 
the cusp profile becomes increasingly flatter. The 
overall mass density drops with an average value 
near a ~ 1.55 since too few heavy mass stars are 
in the cusp to dominate the profile. 

with the same IMF as the other clusters, but did 
not contain an IMBH (dotted lines in lower left 
panel) . The average mass of stars in the core for 
the N — 64K cluster without a black hole rises as 
the cluster evolves towards core collapse (reached 
at T=10.5 Gyrs). Due to the choice of IMF, there 
are only a few black holes present in this model, 
so by the time core collapse has been reached the 
average mass of stars in the core is roughly equal 
to the mass of the most heavy white dwarfs and 
neutron stars. This mass stays nearly constant in 
the post-collapse phase. 

In contrast, in clusters with intermediate-mass 
black holes the average mass of stars in the core 
reaches <m>= 0.6Mq after 2-3 Gyrs and stays 
more or less constant afterwards. There is no de- 
pendence of the average core mass on the particle 
number and no time evolution. This implies that 
our clusters have reached an equilibrium state in 
which heavy stars sink into the cusp due to mass 
segregation and are expelled equally rapidly from 



the cluster centre by a balancing process. The 
most likely driving force for this balancing process 
are close encounters between stars in the cusp. In 
paper I it was shown that close encounters are ef- 
ficient in removing stars from the cusp since the 
average velocities are high, so stars scattered out 
of the cusp leave the cluster completely or are scat- 
tered into the halo and take a long time reaching 
the core again. The average stellar mass in the 
core is similar for star clusters with IMBHs and 
pre-coUapse star clusters without IMBHs, mak- 
ing a distinction between both types of clusters 
through star counts difficult. 

Fig. 7 shows the density distribution of stars 
of different mass groups for the cluster with N — 
131,072 stars. Within the uncertainties, density 
distributions of different mass groups can be fitted 
by power-laws. The exponents a decrease from the 
heavy-mass to the lighter stars since mass segre- 
gation enhances the density of heavy-mass stars 
in the centre. If we exclude the few stellar-mass 
black holes, the most massive stars are massive 
white dwarfs which have masses m < 1.2Mq. 
Fig. 7 shows that their density distribution fol- 
lows a power-law with exponent a — 1.68. For 
the N = 32K and 64K clusters we obtain simi- 
lar exponents of a = 1.80 and a — 1.82 for the 
heavy mass stars. Given the error with which a 
can be determined, all values are probably com- 
patible with an exponent of a — 1.75. Bahcall 
& Wolf (1977) showed that in a two-component 
system, stars of the heavier mass group follow a 
power-law distribution with a — 1.75 in the cusp 
around the black hole. Our simulations are an ex- 
tension of their work to multi-component systems. 
Here again the heaviest stars follow an a = 1.75 
law. Since the mass in the heaviest mass group is 
only a small fraction of the total cusp mass, the 
actual density distribution of the cusp is flatter 
than an a = 1.75 power-law. The slope a for stars 
of average mass m can be approximated by: 



^{m) 



= 0.75 + to/1.1 



(5) 



We find that this density law is a good fit to all 
runs. Tremaine et al. (1994) showed that the slope 
a of the stellar density distribution around a black 
hole has to fulfil the condition a > 1/2 if the ve- 
locity distribution is isotropic, since otherwise the 
average velocity of stars at a given radius r in the 
cusp would become larger than the escape veloc- 



ity at this radius. The density profile of the lowest 
mass stars is indeed above this limit even if the ar- 
gument is not a strict argument for our case since 
Tremaine et al. considered only the escape from 
the potential of the black hole while here the main 
cluster mass is in the form of stars. 

3.4. Tidal disruption of stars 

Fig. 8 shows the relative fraction of disrupted 
stars in the simulations for clusters with parti- 
cle numbers between 16,384 < N < 131,072. 
The stars in our simulations fall roughly into four 
different categories: main-sequence stars, giants, 
white dwarfs and neutron stars and black holes. 
Mostly main-sequence stars are disrupted since 
they are abundant and have relative large radii. 
Giants also contribute significantly despite their 
low numbers. Our simulations show that most gi- 
ants are disrupted within the first Gyr. There are 
no significant differences in the disruption rates 
between individual runs. Taking the mean over all 
simulations, we find that 72%±7% of all disrupted 
stars are main-sequence stars, 19%±3% giants and 
9% ± 2% white dwarfs. During our simulations 
no neutron stars or black holes merged with the 
central IMBH. However, the merging rate of neu- 
tron stars and black holes could be underestimated 
since our simulations did not include gravitational 
radiation, the effects of which will be discussed in 
section 3.6. 

Using the loss-cone theory developed by Frank 
& Rees (1976), one can show with an argumenta- 
tion similar to the one in paper I that the rate at 
which stars are disrupted by a central black hole 
is proportional to 

D^VGl '^ ."°"! ) (6) 



M 



5a — 6 
BH 



where tt is the tidal radius of stars, m their av- 
erage mass and tiq a constant which describes the 
number density of stars in the central cusp around 
the black hole n{r) ~ nor^"". With the help of eq. 
1, we can rewrite this as: 



D = kD^/G 
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In globular clusters, the radii and masses of stars 
differ, so disruption rates have to be calculated 
separately for each stellar type. 
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Fig. 8. — Relative fraction of disrupted stars for 
cluster simulations with different particle num- 
bers. IMBHs in globular clusters preferentially 
disrupt main-sequence stars and giants, with white 
dwarfs accounting for only 10% of all disrupted 
stars. The number of disrupted neutron stars is 
negligible small. 



By comparing the actual number of disruptions 
and the estimated number from eq. 7, we can de- 
termine the coefficient ko- In order to do this, we 
calculated no for all times data was stored from 
the stars inside the cusp, and integrated eq. 7 over 
time to obtain the expected number of disruptions 
for each run and each stellar species. In order to 
calculate no, we assumed a — 1.55 for main se- 
quence stars and giants and a = 1.75 for compact 
remnants. The disruption constants kjj and the 
amount by which different stellar species should 
contribute can be found in Table 2. Most simu- 
lations are compatible with ko — 65, which was 
also found in the single-mass runs of paper I. In 
agreement with the simulations, we expect most 
disrupted stars to be main-sequence stars. Giants 
and white dwarfs should account for about 25% 
of all disruptions and the rate at which neutron 
stars are disrupted is negligible, in agreement with 
the fact that no such disruptions were observed in 
our TV-body runs. The relative fraction of dis- 
rupted giants and white dwarfs agrees very well 
with the results of our runs. Most stars are dis- 
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rupted within the first two Gyrs, when the clusters 
are compact and the densities around the IMBHs 
are high. 

In order to apply our disruption rates to real 
globular clusters, for which the density hq of the 
central cusp is unknown, we have to connect hq 
to the core density n^. We assume that the cusp 
density rises with a = 1.55 and that the cusp goes 
over into a constant density core with density Uc 
at the influence radius of the black hole, given by 
eq. 3. With this, the disruption rate becomes 



D 



1.13 



R^ 



lOOMyrs \R(, 



f Mbh 



m 

aTqJ VioooMq 



105 pC- 



Vr 



10 km/sec 



(8) 



Using this formula and the velocities and cen- 
tral densities of globular clusters given in Pryor 
& Meylan (1993), we can calculate the disruption 
rates D for globular clusters. We assume an aver- 
age stellar mass in the core of m = 0.6Mq and a 
stellar radius i?, — 0.7 Rq. The results are shown 
in Fig. 9. Core-collapse clusters have D values 
of up to several 10"^ yr-i due to their high cen- 
tral densities. However, their morphology rules 
out intermediate-mass black holes, as was shown 
in section 3.2. Among the non-core collapse clus- 
ters, three clusters (NGC 1851, NGC 4147 and 
NGC 6336) have D > 10"^ yr"!, so a 1000 Mq 
IMBH could double its mass within a Hubble time. 
The disruption rate is also high enough that X-ray 
flares from material falling onto the central black 
hole could be observable, depending how long it 
takes for the material to be swallowed by the cen- 
tral black hole. These clusters are, however, more 
concentrated than what would we expect if they 



Table 2: Merger rates of different stellar species. 
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Fig. 9. — Disruption rates D for globular clusters 
from the list of Pryor & Meylan (1993) for which 
central velocity dispersions and densities are avail- 
able. Core-collapse clusters have high disruption 
rates but their morphology shows that they can- 
not contain IMBHs. All clusters which follow the 
mass-radius relation predicted by our runs have 
only small disruption rates, making it unlikely 
that the IMBHs would be visible as X-ray sources. 

contain intermediate-mass black holes. Clusters 
which have sizes that agree with our predicted 
ones have fairly small disruption rates oiD < 10~^ 
yr~^, i.e. less than one star per Gyr. The black 
holes in such clusters would therefore remain in- 
active for most of the time. If we start from an 
50 Mq seed black hole, the disruption rates would 
drop to less than 2.5 • 10^^° yr~^ according to eq. 
8 even for the most favourite clusters. Therefore, 
the growth of an IMBH from a black hole produced 
by normal stellar evolution through the tidal dis- 
ruption of stars is impossible, unless the cluster is 
initially significantly more concentrated than any 
cluster we see today. 

Fig. 10 shows the semi-major axis distribution 
of stars disrupted by the IMBH on their last or- 
bit prior to disruption. Shown is the combined 
distribution for the N = 64K and N = 128K clus- 
ters which are closest to real globular clusters. All 
stars have semi-major axis a which are far larger 
than their tidal radii, similar to the situation found 
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Fig. 10. — Semi-major axis distribution a of stars 
disrupted by the black hole in the iV = 65, 536 and 
N — 131, 072 cluster simulations. Most stars move 
on highly eccentric orbits with a » rT when 
they are disrupted by the black hole. Disrupted 
white dwarfs move on orbits with semi-major axes 
that are about a factor of 10 smaller than main- 
sequence stars. 

for single-mass clusters in paper I. The most im- 
portant mechanism for the disruption of stars is 
therefore again drift in angular momentum space. 
Material from disrupted stars will first move in a 
highly eccentric ring around the IMBH. The ring 
shrinks due to viscous heating until the gas from 
the disrupted star forms an accretion disc around 
the IMBH. Since there are other stars moving in- 
side the initial orbit of the disrupted star which 
can scatter away material or swallow it, the frac- 
tion of material that is finally swallowed by the 
IMBH is rather uncertain. 

Due to their smaller radii, disrupted white 
dwarfs also come from smaller distances. Ac- 
cording to the theory developed in paper I for 
single-mass clusters, their average semi-major axis 
should be smaller by about {Rms/Rwd)"^^^ ~ 10, 
which is in good agreement with the results in our 
simulations. Similarly, giants come from larger 
distances on average. 
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Fig. 11. — Number of bound black holes with 
masses Mbh > 5Mq in the runs with a high up- 
per mass limit. The number decreases due to close 
encounters between the black holes in the central 
cusp around the IMBH. After several Gyrs, only 
one heavy mass black hole remains in most cases. 
This is the object most deeply bound to the cen- 
tral IMBH and has absorbed the energy from the 
encounters which ejected the other black holes. 



3.5. Clusters with high-mass black holes 

We now discuss the evolution of clusters which 
started with a mass function extending up to 100 
Mq. These clusters contained a significant num- 
ber of stellar mass black holes. The highest mass 
black holes formed had about 45 Mq and we as- 
sumed a 100% black hole retention rate, so the 
true fraction of black holes in globular clusters will 
probably be somewhere between the situation in 
these runs and the previous ones. Cluster radii 
were chosen such that the relaxation time was the 
same in all runs and equal to the relaxation time 
of a dense globular cluster with Mc = 10^ Mq and 
rh ^ I pc. 

Massive black holes are initially produced 
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throughout the clusters and sink towards the cen- 
tres due to mass segregation. Fig. 11 shows the 
number of black holes bound to the clusters as a 
function of time. The number of heavy stellar- 
mass black holes drops since close encounters 
between the black holes in the cusp around the 
IMBH remove them from the cluster. Since the 
velocity in the cusp is relatively high, such en- 
counters mostly lead to the ejection of one of the 
black holes. For the N = 64K and N = 128K 
clusters, only one massive black hole remains in 
the cluster after several Gyrs have passed. For 
the N = 32K cluster 3 remained, however two 
of them were in orbits far away from the cluster 
centre and would have been lost from the cluster 
if the cluster would have been surrounded by a 
galactic tidal field. In all three clusters, the black 
hole which remains is the star closest bound to the 
central IMBH (see Fig. 14) and is among the most 
heavy black holes produced. This resembles the 
situation in our runs with the lower upper mass 
limit. As the relaxation time in our runs is similar 
to that of dense, massive globular clusters, we ex- 
pect that clusters with IMBHs also contain only 
few other massive black holes, and the star closest 
bound to the IMBH should be another black hole. 
A look at the lagrangian radii (Fig. 12) shows 
that the overall expansion is slightly different from 
the previous case. Clusters expand more rapidly in 
the beginning when they still contain many heavy 
mass black holes. These are very effective in scat- 
tering low-mass stars to less bound orbits while 
sinking into the cluster centre. The innermost 
radii of the cluster with N = 128K stars and a high 
upper mass limit decrease slightly after T = 7000 
Myrs. This is due to the loss of the 2nd nearest 
star to the IMBH through a close encounter with 
the innermost star (see Fig. 14). The 2nd near- 
est star was an efficient heat source since it was 
a 15 Mq black hole which moved in a relatively 
wide orbit around the IMBH, bringing it into fre- 
quent encounters with field stars. After this BH 
is lost, the innermost radii shrink to adjust them- 
selves and the energy generation rate in the centre 
to the size of the half-mass radius. The final radii 
of clusters with heavy mass black holes are within 
10% of the radii of runs with the same N but a 
lower upper mass limit (see Table 1). The ini- 
tial radii were different for N > 64K stars, but it 
was shown in section 3.1 that the initial radii do 
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Fig. 12. - Evolution of lagrangian radii for two 
N = 128K clusters, one having an IMF extend- 
ing up to ?Ti = IOOMq and many high-mass black 
holes and one with a lower high-mass cut-off of 30 
Mq. The cluster with many black holes expands 
stronger initially due to the efficient heating of the 
black holes. Nevertheless, the final radii are al- 
most the same for both clusters. 



not influence the final radii much. We therefore 
conclude that the initial mass function does not 
significantly affect the final radius either. 

Fig. 13 depicts the density distribution of the 
A^ — 128K cluster after 12 Gyrs. As in the case of 
a cluster starting with only few black holes, high- 
mass stars follow a steeper density distribution 
and are enriched in the cluster centre. The overall 
slopes for stars of the same mass are very similar 
to the low upper mass limit case since most high- 
mass black holes have been lost from the cluster 
by this time, so the overall mass function of stars 
is nearly identical. Again, when viewed in projec- 
tion, this cluster would appear as a cluster with a 
constant density core. The average mass of stars 
in the core is the same as in the previous case, 
<m>= O.QMq, independent of N. 

The merger rates of different stars for clusters 
with a high upper mass limit can be found in Ta- 
ble 2. Compared to the lower mass limit case, 
the fraction of neutron star and black hole dis- 
ruptions is increased by a factor of 100 due to 
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Fig. 13. — Density of different mass groups as 
function of radius for the A^ = 128K star cluster 
with a high upper mass cut-off at T=12 Gyrs. The 
power-law slopes are similar to that of the cluster 
with a low upper mass limit. High-mass stars are 
again enriched in the cluster centre. 

their higher density around the IMBH in the ini- 
tial phases. The overall fraction is however still 
negligible. The merger rate for giants decreased 
since the high-mass stellar mass black holes pre- 
vent giants from accumulating near the IMBH in 
the initial phases. The effect is less visible for the 
N = 128K clusters, in which the disruption rate is 
already quite similar to the run with the low upper 
mass limit. The constant for the overall disruption 
rate ko stayed more or less the same, so our pre- 
vious conclusions still hold if we change the mass 
function. This is despite the fact that the IMBH 
is forced to move with a larger amplitude since the 
star closest bound to the IMBH is now 20 times 
more massive. The movement of the central black 
hole has therefore not much effect on the merging 
rate. 

3.6. Gravitational radiation 

Figs. 14 and 15 depict the semi-major axis of 
stars which are most deeply bound to the IMBH 
for the two clusters with N = 128K stars. The en- 
ergy of the deepest bound star decreases quickly 
in the beginning when it still has many interac- 



tions with passing stars. When the semi-major 
axis becomes significantly smaller than that of 
other stars, interactions become rare and the en- 
ergy change slows down considerably. In both 
clusters, the innermost star is among the heavi- 
est stars in the cluster, and would be a black hole 
with several IOMq for a globular cluster with a 
reasonable IMF. The innermost star will therefore 
not transfer mass onto the IMBH. All other stars 
have semi-major axis oi R > lO^i?0 which is too 
far for mass transfer, even if some stars will move 
on strongly radial orbits. 

The time for two black holes in orbit around 
each other to merge due to emission of gravita- 
tional radiation is equal to (Evans et al. 1987): 
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where a is the semi-major axis of the orbit, c the 
speed of light, m\ and m^ are the masses of the 
two black holes, M the combined mass and F is a 
function of the orbital eccentricity e and is given 
by (Peters 1964): 
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The eccentricity of the orbit of the innermost star 
fluctuates rapidly as long as the star undergoes 
many close encounters with passing stars, but be- 
comes fairly stable as soon as the star has detached 
itself from the other cluster stars. Most of the 
time the innermost stars moves on an orbit with 
moderate eccentricity and we can assume e = 0.5. 
With this value, the radius at which a black hole 
of 20Mq merges with a 1000 Mq IMBH within a 
Hubble time is a = 562i?Q. This radius is marked 
by a dashed line in Figs. 14 and 15. The semi- 
major axis for the innermost stars in our clusters 
are roughly a factor of 6 higher at the end of the 
simulation, so the merging timescale is around a 
thousand Hubble times, too long to have any no- 
ticeable effect on the orbits. 

For a circular orbit, the frequency / and ampli- 
tude h of gravitational waves emitted by two black 
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holes in orbit around each other at a distance R 
are given by (Douglas & Braginsky 1979): 
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Fig. 16 shows the amplitude and frequency of grav- 
itational radiation emitted from the N — 128K 
clusters which started with a high upper mass 
limit. Each time the data were stored, we cal- 
culated both values, assuming that the cluster is 
at a distance of i? = 8 kpc. The innermost black 
hole is in too wide an orbit to be detectable. 

How will these results change for globular clus- 
ters ? In clusters with higher particle numbers, the 
distance of the innermost stars will be different. In 
paper I it was shown that the energy generation 
rate in the cusp is proportional to: 
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for an a = 1.75 cusp. Such a profile should be 
established very close to the IMBH in higher-A^ 
models. The rate at which energy can be trans- 
fered outward at the half-mass radius is given by: 
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if the weak N dependence in the Coulomb loga- 
rithm can be neglected. Since for a cluster evolv- 
ing slowly along a sequence of equilibrium models 
energy generation must balance energy transport, 
a condition for tiq can be obtained from both equa- 
tions. The distances of the innermost stars from 
the IMBH in an a = 1.75 cusp follow a relation 
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Fig. 14. — Semi-major axis of the three stars deep- 
est bound to the IMBH as function of time for 
the cluster with TV = 128K stars and a high up- 
per mass limit. The object closest bound to the 
IMBH is almost always another black hole which 
is among the heaviest objects in the cluster. The 
other stars are too far away from the BH to un- 
dergo mass transfer. 
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Tj ~ r7g . We therefore obtain for the distance 
of the innermost star from the IMBH: 
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Using the relation found for rh in section 3.1, rh ~ 

— 1/3 

M(-, , we obtain for the dependence of r^ on the 
cluster mass Mc'- 



M, 



-0.93 



(16) 



i.e. the distances of the innermost stars decrease 
nearly linear with the particle number. Assuming 
an a = 1.55 cusp gives nearly the same scaling 
law. The data in our A'^-body runs is consistent 
with this relation. Eq. 16 predicts that even in 
the most massive globular clusters, the distance of 
the innermost main-sequence stars from the IMBH 
is larger than 10'^ Rq, too large for stable mass 
transfer. Dynamical evolution alone is therefore 



15 



10' 

m 

K. 104 

1000 
10' 

1 10= 

ID 
n 

« lO' b- 



^f\i^^\{fH^'f'^^^m^^ 



Main sequence star 
Red giant 



Compact remnant 



H ! \ H 



iM#f#1|iV^ 



H ^ 1 \ 




5000 

T [Myrs] 

Fig. 15. — Same as Fig. 14 for the case of a low 
upper mass limit. In this case, the star closest 
bound to the IMBH is a heavy-mass white dwarf 
or neutron star. 



not strong enough to form stable X-ray binaries 
involving the IMBH. 

As Hopman et al. (2004) have shown, it might 
be possible for an IMBH to capture a passing star 
through tidal heating. In this case the star could 
end up in a circular orbit with small enough radius 
so that mass transfer onto the IMBH is possible. 
Encounters with passing stars will however scat- 
ter the star either out of the cusp or onto a highly 
eccentric orbit where it is disrupted by the IMBH. 
The present paper does not include orbital changes 
due to the tidal heating of stars, so detailed simu- 
lations have to be done to study the tidal capture 
rate of an IMBH and the further orbital evolution 
of the captured stars. Apart from stars, the IMBH 
could also accrete cluster gas lost from post-main 
sequence stars through stellar winds. The result- 
ing X-ray flux depends on the gas fraction in the 
cluster and the details of the accretion mechanism, 
but could be bright enough to be observable under 
favourable conditions (Ho, Terashima & Okajima 
2003). 

The situation looks more promising for gravi- 
tational radiation. A factor of 10 increase in the 
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Fig. 16. — Frequency and amplitude of gravita- 
tional radiation emitted from the innermost BH 
binary for the N = 128K star cluster which starts 
with a high upper mass limit. Although the bi- 
nary hardens as a result of encounters with pass- 
ing stars, the orbital separation is still too large, 
and hence the frequency too small, for the binary 
to become detectable. 



particle number from our largest N — 128K star 
runs would already be enough that a tight IMBH- 
BH binary forms which merges within a Hubble 
time. If the semi-major axis of a 10 Mq-IOOO 
M© BH binary in a galactic globular cluster is 
less than about 50 Rq, it would become visible. 
By this time, the time remaining to final merg- 
ing has dropped to 10® yrs. Making the conserva- 
tive assumption that each globular cluster with an 
IMBH goes only once through such a phase, and 
assuming that 10% of all galactic globular clusters 
contain IMBHs, chances are around 10% that any 
galactic globular cluster is presently emitting de- 
tectable amounts of gravitational radiation. If a 
GW source has a signal-to-noise level larger than 
2, LISA will have an angular resolution of a few 
degrees, making it possible to identify the globular 
cluster containing the IMBH. 

The detection is even more likely if we con- 
sider extragalactic globular clusters. When the 
semi-major axis has dropped to a = IRq, a 10 
M0-1OOO Mq BH binary would be bright enough 
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to be visible to distances of i? « 1 Gpc. Based 
on the 2dF galaxy redshift survey, Norberg et al. 
(2002) estimated a luminosity density of pL ~ 
(1.8 ± 0.17) • lO^L© Mpc~3 at z = 0, which cor- 
responds to roughly 0.1 Milky Way sized galaxies 
every Mpc^^. Assuming that each of them con- 
tains 100 globular clusters gives 5 • 10^" globular 
clusters inside 1 Gpc. Assuming again that 10% of 
all globular clusters contain an IMBH and that the 
merging rate of black holes with the central IMBH 
is constant over time, gives 5 events that would 
be visible with LISA at any given time. Globular 
clusters containing IMBHs will therefore be an im- 
portant source of GW emission from star clusters, 
in addition to double compact stars (Benacquista 
et al. 2001). 

4. Conclusions 

We have performed a set of large iV-body sim- 
ulations of multi-mass star clusters containing 
intermediate-mass black holes. Our simulations 
include a realistic mass spectrum of cluster stars, 
mass-loss due to stellar evolution, two-body relax- 
ation and tidal disruption of stars by the central 
black hole. These simulations are the first fully 
self-consistent simulations of realistic star clusters 
with IMBHs. 

Our results can be summarised as follows: A 
density cusp forms around the central black hole 
with a density profile p ^ 7.-1-55 jj-^ three dimen- 
sions. For low-mass IMBHs with a mass less than 
a few percent of the cluster mass, the cusp extends 
out to a radius where the velocity dispersion in the 
core becomes comparable to the circular velocity 
of stars around the black hole. In this case, the 
stars in the cusp contain only a fraction of the 
mass of the central black hole, which makes the 
direct detection of the cusp difficult for IMBHs of 
Mbh < 1OOOM0. Globular clusters with IMBHs 
following the relation found by Gebhardt et al. 
(2000) for galactic bulges belong to this category. 
Only more massive black holes create a power-law 
cusp profile throughout the cluster core that would 
be directly visible. 

When viewed in projection, the luminosity 
profiles of clusters with massive black holes dis- 
play a constant density core. The presence of 
intermediate-mass black holes in galactic core- 
collapse clusters like Ml 5 is therefore ruled out 



unless these clusters have a stellar mass distri- 
bution very different from our clusters. As was 
shown in Baumgardt et al. (2003a), a more nat- 
ural explanation for mass-to-light ratios that in- 
crease toward the centre in such clusters is a dense 
concentration of neutron stars, white dwarfs and 
stellar mass black holes. The amount of mass 
segregation in a cluster with an IMBH is also 
smaller compared to a post-core collapse cluster, 
the average mass of stars in the centre being about 
m = 0.6 Mq. Clusters with IMBHs therefore re- 
semble star clusters which are in the pre-coUapse 
phase also in terms of the amount of mass segre- 
gation. 

All clusters with intermediate-mass black holes 
expand due to close encounters of stars in the cusp 
around the central black hole. We find that the 
values of the half-mass radii reached depend on 
the mass of the central black hole and the number 
of cluster stars, but are nearly independent of the 
initial cluster radius and density profile. Portegies 
Zwart et al. (2004) have shown that central den- 
sities of more than lO^M^g/pc'^ are necessary to 
form an IMBH through runaway merging of mas- 
sive main-sequence stars. Similarly high densities 
are necessary to form an IMBH through the merg- 
ing of stellar mass black holes through gravita- 
tional radiation (Mouri & Taniguchi 2002). Such 
densities are among the highest found in globular 
clusters. Our simulations show that even if clus- 
ters with intermediate-mass black holes start with 
very high densities, the subsequent cluster expan- 
sion is sufficient to put them among the least con- 
centrated clusters after a Hubble time. Low-mass 
clusters surrounded by a strong tidal field will dis- 
solve due to the cluster expansion, releasing their 
IMBHs. For clusters close enough to a galactic 
centre, these IMBHs could then spiral into the 
centre and merge through the emission of gravi- 
tational radiation. If enough IMBHs are formed, 
this process might provide the seed black holes for 
the supcrmassive black holes observed in galactic 
centres (Ebisuzaki et al. 2001). 

IMBHs in star clusters disrupt mainly main- 
sequence stars and giants, with white dwarfs ac- 
counting for only 10% of all disruptions. In young 
star clusters, the tidal disruption rate of giants 
is similar to that of main-sequence stars. During 
our simulations, no neutron stars were disrupted, 
so intermediate-mass black holes in star clusters 
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do not emit gamma rays or gravitational radiation 
from such events. Most stars that were disrupted 
moved around the central black hole on highly ec- 
centric orbits with large semi-major axis. Even 
if 100% of the mass from disrupted stars is being 
accreted onto the central black hole, tidal disrup- 
tions of stars are too rare to form an intermediate- 
mass black hole out of a Mbh ^ 50Mq progen- 
itor, except the initial cluster was significantly 
more concentrated than present day globular clus- 
ters. In such cases, the density would however also 
be high enough that run-away merging of heavy 
mass main-sequence stars would lead directly to 
the formation of an IMBH. The largest disruption 
rates for globular clusters which do not have core- 
collapse profiles are about 10~^/yr for clusters 
with half-mass radii significantly smaller than pre- 
dicted by our runs. Clusters with half-mass radii 
in agreement with our simulations have two or- 
ders of magnitudes smaller disruption rates. The 
black holes in such clusters are therefore inactive 
for most of the time. 

The detection of a IOOOMq IMBH in a globular 
cluster through the measurement of radial veloci- 
ties or proper motions of cluster stars requires the 
observation of about 20 stars in the central cusp. 
For globular clusters that are close enough, the 
central cusp extends to distances of several arc- 
sec, so the detection of an intermediate-mass black 
hole should be possible by either radial velocity or 
proper motion studies with HST. 

Black holes with masses of SM© or higher are 
strongly depleted in star clusters with intermediate- 
mass black holes since they sink into the cen- 
tre through dynamical friction and then remove 
each other by close encounters in the central cusp 
around the IMBH (Kulkarni et al. 1993). Based 
on our simulations, we expect that only one high- 
mass black hole remains in the cluster. This black 
hole is among the heaviest black holes formed and 
ends up as the object most tightly bound to the 
IMBH. In our runs, the distance of the innermost 
black hole to the IMBH was never small enough 
that the frequency of gravitational radiation was 
in the range observable for e.g. LISA. For higher 
particle numbers or for clusters that start off more 
concentrated, it seems likely that a tight enough 
BH-IMBH binary is formed dynamically. If 10% 
of all globular clusters contain IMBHs, and each 
IMBH merges with a stellar-mass black hole at 



least once within a Hubble time, the chance that 
any galactic globular cluster currently has a tight 
enough BH-IMBH binary detectable for LISA is 
10%. Within a radius R= \ Gpc, about 5 globu- 
lar clusters would harbour bright enough sources 
for LISA at any one time. 
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